Search for: All records

ABSTRACT Theoretical arguments and observations suggest that in massive haloes ($$>10^{12}\, {\rm M}_\odot$$), the circumgalactic medium (CGM) is dominated by a ‘hot’ phase with gas temperature near the virial temperature ($$T\approx T_{\rm vir}$$) and a quasi-hydrostatic pressure profile. Lower-mass haloes are however unlikely to be filled with a similar quasi-static hot phase, due to rapid radiative cooling. Using the FIRE (Feedback In Realistic Environment) cosmological zoom simulations, we demonstrate that the hot phase is indeed subdominant at inner radii ($$\lesssim 0.3 R_{\rm vir}$$) of $$\lesssim 10^{12}\, {\rm M}_\odot$$ haloes, and the inner CGM is instead filled with $$T\ll T_{\rm vir}$$ gas originating in outflows and inflows, with a turbulent velocity comparable to the halo virial velocity. The turbulent velocity thus exceeds the mass-weighted sound speed in the inner CGM, and the turbulence is supersonic. UV absorption features from such CGM trace the wide lognormal density distributions of the predominantly cool and turbulent volume-filling phase, in contrast with tracing localized cool ‘clouds’ embedded in a hot medium. We predict equivalent widths of $$W_\lambda \sim 2\lambda v_{\rm c}/c\sim 1$$Å for a broad range of strong UV and EUV transitions (Mg ii, C ii, C iv, Si ii–iv, O iii–v) in sightlines through inner CGM dominated by turbulent pressure of $$\lesssim L^\star$$ galaxies at redshifts $$0\le z\lesssim 2$$, where $$\lambda$$ is the transition wavelength, $$v_{\rm c}$$ is the circular velocity, and c is the speed of light. Comparison of our predictions with observational constraints suggests that star forming $$\lesssim$$ $$L^\star$$ and dwarf galaxies are generally dominated by turbulent pressure in their inner CGM, rather than by thermal pressure. The inner CGM surrounding these galaxies is thus qualitatively distinct from that around quenched galaxies and massive discs such as the Milky-Way and M31, in which thermal pressure likely dominates. 

ABSTRACT The cold ($$\sim 10^{4}\, {\rm K}$$) component of the circumgalactic medium (CGM) accounts for a significant fraction of all galactic baryons. However, using current galaxy-scale simulations to determine the origin and evolution of cold CGM gas poses a significant challenge, since it is computationally infeasible to directly simulate a galactic halo alongside the sub-pc scales that are crucial for understanding the interactions between cold CGM gas and the surrounding ‘hot’ medium. In this work, we introduce a new approach: the Cold Gas Subgrid Model (CGSM), which models unresolved cold gas as a second fluid in addition to the standard ‘normal’ gas fluid. The CGSM tracks the total mass density and bulk momentum of unresolved cold gas, deriving the properties of its unresolved cloudlets from the resolved gas phase. The interactions between the subgrid cold fluid and the resolved fluid are modelled by prescriptions from high-resolution simulations of ‘cloud crushing’ and thermal instability. Through a series of idealized tests, we demonstrate the CGSM’s ability to overcome the resolution limitations of traditional hydrodynamics simulations, successfully capturing the correct cold gas mass, its spatial distribution, and the time-scales for cloud destruction and growth. We discuss the implications of using this model in cosmological simulations to more accurately represent the microphysics that govern the galactic baryon cycle. 

Abstract We presentCloudFlex, an open-source tool for predicting absorption-line signatures of cool gas in galaxy halos with small-scale structure. Motivated by analyses of ∼10⁴K material in hydrodynamical simulations of turbulent, multiphase media, we model cool gas structures as complexes of cloudlets sampled from a power-law distribution of mass $\propto m_{\mathrm{cl}}^{-\alpha }$with velocities drawn from a turbulent velocity field. The user may specifyα, the lower limit of the cloudlet mass distribution ( $m_{\mathrm{cl},\min }$), and several other parameters that set the mass, size, and velocity distribution of the complex. This permits investigation of the relation between these parameters and absorption-line observables. As a proof-of-concept, we calculate the Mgiiλ2796 absorption induced by the cloudlets in background quasi-stellar object (QSO) spectra. We demonstrate that, at fixed metallicity, the covering fraction of sight lines with equivalent widthsW₂₇₉₆< 0.3 Å increases significantly with decreasing $m_{\mathrm{cl},\min }$, cloudlet number density (n_cl), and complex size. We then use this framework to predict the halo-scaleW₂₇₉₆distribution around ∼L^*galaxies. We show that the observed incidences ofW₂₇₉₆> 0.3 Å sight lines with impact parameters 10 kpc <R_⊥< 50 kpc in projected QSO–galaxy studies are consistent with our model over much of parameter space. However, they are underpredicted by models with $m_{\mathrm{cl},\min }\ge 100M_{\odot }$andn_cl≥ 0.03 cm⁻³, in keeping with a picture in which the inner cool circumgalactic medium (CGM) is dominated by numerous low-mass cloudlets (m_cl≲ 100M_⊙) with a volume filling factor ≲1%. When used to model absorption-line data sets built from multi-sight line and/or spatially extended background probes,CloudFlexenables detailed constraints on the size and velocity distributions of structures comprising the photoionized CGM. 

Abstract The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the circumgalactic medium (CGM) is expected to result in a bow shock that leads the LMC’s gaseous disk. In this letter, we use hydrodynamic simulations of the LMC’s recent infall to predict the extent of this shock and its effect on the Milky Way’s (MW) CGM. The simulations clearly predict the existence of an asymmetric shock with a present-day standoff radius of ∼6.7 kpc and a transverse diameter of ∼30 kpc. Over the past 500 Myr, ∼8% of the MW’s CGM in the southern hemisphere should have interacted with the shock front. This interaction may have had the effect of smoothing over inhomogeneities and increasing mixing in the MW CGM. We find observational evidence of the existence of the bow shock in recent Hαmaps of the LMC, providing a potential explanation for the envelope of ionized gas surrounding the LMC. Furthermore, the interaction of the bow shock with the MW CGM may also explain the observations of ionized gas surrounding the Magellanic Stream. Using recent orbital histories of MW satellites, we find that many satellites have likely interacted with the LMC shock. Additionally, the dwarf galaxy Ret2 is currently sitting inside the shock, which may impact the interpretation of the reported gamma-ray excess in Ret2. This work highlights how bow shocks associated with infalling satellites are an underexplored yet potentially very important dynamical mixing process in the circumgalactic and intracluster media. 

ABSTRACT Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological fire-2 CR-magnetohydrodynamic simulations of L* galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics-driven hysteresis. SC models show a higher tendency to undergo ‘ejective’ feedback events due to a runaway buildup of CR pressure in dense gas due to the behaviour of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, interstellar medium gas, and B properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L* galaxies, and indicating the prospect for non-thermal radio continuum observations to constrain CR transport physics. 

Abstract Synchrotron emission is one of few observable tracers of galactic magnetic fields (B) and cosmic rays (CRs). Much of our understanding of B in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however it remains unclear how well these assumptions hold, and what B these estimates physically represent. Using FIRE simulations which self consistently evolve CR proton, electron, and positron spectra from MeV to TeV energies, we present the first synthetic synchrotron emission predictions from simulated L* galaxies with ”live” spectrally-resolved CR-MHD. We find that synchrotron emission can be dominated by relatively cool and dense gas, resulting in equipartition estimates of B with fiducial assumptions underestimating the ”true” B in the gas that contributes the most emission by factors of 2-3 due to small volume filling factors. Motivated by our results, we present an analytic framework that expands upon equipartition models for estimating B in a multi-phase medium. Comparing our spectrally-resolved synchrotron predictions to simpler spectral assumptions used in galaxy simulations with CRs, we find that spectral evolution can be crucial for accurate synchrotron calculations towards galactic centers, where loss terms are large. 

Abstract We make an in-depth analysis of different active galactic nuclei (AGN) jet models’ signatures, inducing quiescence in galaxies with a halo mass of 10¹²M_⊙. Three jet models, including cosmic-ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. Each of our simulations is idealized isolated galaxy simulations with AGN jet powers that are constant in time and generated using GIZMO and with FIRE stellar feedback. We examine the distribution of Mgii, Ovi, and Oviiiions, alongside gas temperature and density profiles. Low-energy ions, like Mgii, concentrate in the interstellar medium (ISM), while higher energy ions, e.g., Oviii, prevail at the AGN jet cocoon’s edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic-ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic-ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in Oviand Oviiiis demonstrated in the ISM-to-circumgalactic medium (CGM) transition (between 10 and 30 kpc), growing smoothly toward the CGM at greater distances. 

ABSTRACT Recent theoretical studies predict that the circumgalactic medium (CGM) around low-redshift, ∼L* galaxies could have substantial non-thermal pressure support in the form of cosmic rays. However, these predictions are sensitive to the specific model of cosmic ray transport employed, which is theoretically and observationally underconstrained. In this work, we propose a novel observational constraint for calculating the lower limit of the radially averaged, effective cosmic ray transport rate, $${\kappa _{\rm eff}^{\rm min}}$$. Under a wide range of assumptions (so long as cosmic rays do not lose a significant fraction of their energy in the galactic disc, regardless of whether the cosmic ray pressure is important or not in the CGM), we demonstrate a well-defined relationship between $${\kappa _{\rm eff}^{\rm min}}$$ and three observable galaxy properties: the total hydrogen column density, the average star formation rate, and the gas circular velocity. We use a suite of Feedback in Realistic Environments 2 galaxy simulations with a variety of cosmic ray transport physics to demonstrate that our analytical model of $${\kappa _{\rm eff}^{\rm min}}$$ is a robust lower limit of the true cosmic ray transport rate. We then apply our new model to calculate $${\kappa _{\rm eff}^{\rm min}}$$ for galaxies in the COS-Halos sample, and confirm this already reveals strong evidence for an effective transport rate that rises rapidly away from the interstellar medium to values $${\kappa _{\rm eff}^{\rm min}}\gtrsim 10^{30\!-\!31}\, {\rm cm}^2\, {\rm s}^{-1}$$ (corresponding to anisotropic streaming velocities of $$v^{\rm stream}_{\rm eff} \gtrsim 1000\, {\rm km}\, {\rm s}^{-1}$$) in the diffuse CGM, at impact parameters larger than 50–100 kpc. We discuss how future observations can provide qualitatively new constraints in our understanding of cosmic rays in the CGM and intergalactic medium. 

ABSTRACT The physics of magnetic fields (B) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyse the magnetic fields in a set of high-resolution, magnetohydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionized interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM |B| fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circumgalactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large (> 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their haloes. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to neutral hydrogen number densities < 300 cm−3, in contrast to observationally derived models, but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10–200 kpc and better constrain theoretical predictions. 

Abstract We present the KODIAQ-Z survey aimed to characterize the cool, photoionized gas at 2.2 ≲z≲ 3.6 in 202 Hi-selected absorbers with 14.6 ≤ $\log N_{\mathrm{H}\mathrm{I}}$< 20 that probe the interface between galaxies and the intergalactic medium (IGM). We find that gas with $14.6\le \log N_{\mathrm{H}\mathrm{I}}<20$at 2.2 ≲z≲ 3.6 can be metal-rich (−1.6 ≲ [X/H] ≲ − 0.2) as seen in damped Lyαabsorbers (DLAs); it can also be very metal-poor ([X/H] < − 2.4) or even pristine ([X/H] < − 3.8), which is not observed in DLAs but is common in the IGM. For $16<\log N_{\mathrm{H}\mathrm{I}}<20$absorbers, the frequency of pristine absorbers is about 1%–10%, while for $14.6\le \log N_{\mathrm{H}\mathrm{I}}\le 16$absorbers it is 10%–20%, similar to the diffuse IGM. Supersolar gas is extremely rare (<1%) at these redshifts. The factor of several thousand spread from the lowest to highest metallicities and large metallicity variations (a factor of a few to >100) between absorbers separated by less than Δv< 500 km s⁻¹imply that the metals are poorly mixed in $14.6\le \log N_{\mathrm{H}\mathrm{I}}<20$gas. We show that these photoionized absorbers contribute to about 14% of the cosmic baryons and 45% of the cosmic metals at 2.2 ≲z≲ 3.6. We find that the mean metallicity increases withN_Hi, consistent with what is found inz< 1 gas. The metallicity of gas in this column density regime has increased by a factor ∼8 from 2.2 ≲z≲ 3.6 toz< 1, but the contribution of the $14.6\le \log N_{\mathrm{H}\mathrm{I}}<19$absorbers to the total metal budget of the universe atz< 1 is a quarter of that at 2.2 ≲z≲ 3.6. We show that FOGGIE cosmological zoom-in simulations have a similar evolution of [X/H] withN_Hi, which is not observed in lower-resolution simulations. In these simulations, very metal-poor absorbers with [X/H] < − 2.4 atz∼ 2–3 are tracers of inflows, while higher-metallicity absorbers are a mixture of inflows and outflows.